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Abstract 

We present a dynamical study on the nonlinear conduction behaviour in the commensurate 
charge-density-wave phase of the quasi-one-dimensional conductor DCNQBCu below 75 K. We 
can accurately simulate magnitude and time-dependence of the measured conductivity in response 
to large voltage pulses by accounting for the energy exchange between the phononic and electronic 
subsystems by means of an electrothermal model. Our simulations reveal a distinct non-equilibrium 
population of optical phonon states with an average energy of Eph = 19 meV being half the ac¬ 
tivation energy of about AEa = 39 meV observed in DC resistivity measurements. By inelastic 
scattering, this hot optical phonon bath generates additional charge-carrying excitations thus pro¬ 
viding a multiplication effect while energy transferred to the acoustic phonons is dissipated out of 
the system via heat conduction. Therefore, in high electric helds a preferred interaction of charge¬ 
carrying excitations with optical phonons compared to acoustic phonon modes is considered to be 
responsible for the nonlinear conduction effects observed in DCNQFCu. 

PACS numbers: 71.30.-^h,71.45.Lr,72.10.Di,72.20.Ht,63.20.kk 
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I. INTRODUCTION 


Organic semiconductors have been widely investigated in recent decades with respect to 
their potential application in electronic devices. Often their electron and hole mobilities 
are limited due to an efficient interaction of charge carriers with phonons the latter being 
lower in energy compared to conventional inorganic semiconductors because of the weak 
intermolecular binding forces. More exotic semiconducting ground states, such as Peierls or 
Mott insulators, occur in low-dimensional organic conductors often leading to nonlinear con¬ 
duction phenomena!^^. These nonlinearities seem to be distinctive for materials, otherwise 
distinguished by quite different ground states, such as spin-Peierls systems like K-TCNQP 
or charge-ordered states as in a-(BEDT-TTF) and therefore aim for a more general 
description. Mori et al.^ were the hrst to recognize the connection between the steepness 
of the resistance curves and the electric threshold helds above which nonlinear conduction 
occurred. They proposed an electrothermal model to describe the current-volt age charac¬ 
teristics explicitly considering electrical (Joule) heating of distinct electronic or phononic 
degrees of freedom in the material. In this model, the effective specihc heat is the crucial 
parameter hinting at the mechanism that governs the nonlinear transport effects. This quan¬ 
tity is usually interpreted as electronic specihc heat caused by the heating of the electronic 
system at concurrent inefficient energy transfer between electrons and phonons in contrast 
to conventional organic semiconductors. 

To study nonlinear conduction and the validity of the electrothermal model suggested by 
Mori et al.El, we chose Copper radical anion salts of the Dimethyl-Dicyanoquinonediimine 
(DCNQI) molecule, hrst synthesized and studied in the group of Hueni^ as an archetypi¬ 
cal low-dimensional molecular metal. Cu and DCNQI form a quasi-one-dimensional radical 
anion salt with separate stacking structure staying metallic down to lowest temperatures. 
In (DCNQI-d6)2Cu the methyl groups are deuterated causing a three-fold commensurate 
lattice distortion along the face-to-face stacked molecular columns and charge ordering of 
same periodicity at the Cu sites below around 80 K. This phase transition from the metallic 
to the commensurate charge-density-wave (CDW) ground state is of hrst order resulting 
in a step-like increase of the conduct ivit}^. Nonlinear conduction in the CDW state of 
DCNQI2CU radical anion salts has been observed beford^^^^. By means of dielectric response 
measurements the excitation of soliton-like charged domain wall pairs has been identihed as 


2 



relaxation channel indicating their importance for the transport properties in the semicon¬ 
ducting phase of DCNQl2Ciin'. 

Furthermore, the ground state of deuterated DCNQI2CU is very sensitive to the change in 
energy when switching from the free electron gas in the metallic phase to localized electrons 
in the insulating phas^^, i.e. a gain in electronic energy or delocalization of electrons by 
an applied electric held should be more likely to induce an insulator-metal transition com¬ 
pared to other systems. Finally, the availability of high-quality literature dataP^ on the 
temperature-dependent contributions of the electronic system as well as of the acoustic and 
optical phonons to the specihc heat allows for an identihcation of the leading microscopic 
mechanism causing nonlinear conduction. 

In this paper we will show that nonlinear conduction in DCNQI2CU can be well described 
by means of the above-mentioned electrothermal model. For this purpose we carried out 
detailed transient analyses of the conduction behavior allowing for a comprehensive deter¬ 
mination of the effective specihc heat. The magnitude of this quantity and its temperature- 
dependence rendered it possible to identify its major contribution from optical phonons with 
an energy of about half the activation energy deduced from the temperature-dependent low- 
held conductivity. Based on these observations we develop an important extension of the 
established electrothermal model by including the distinct electrical heating of the optical 
phonon subsystem requiring an efficient interaction of the charge-carrying excitations with 
optical phonons similar to conventional organic semiconductor^. 


II. THEORY 

In order to account for the universal nonlinear conduction behavior in a variety of organic 
charge-transfer materials with diherent ground states, Mori et ah proposed an electrothermal 
modeM 


nC— = P + V ■ {kVT} ( 1 ) 

where nC denotes the ehective heat capacity per volume, P is the supplied electrical 
energy, k is the thermal conductivity and T is the temperature of the system. The supplied 
electrical energy reads P = a{T)E‘^, where cr(T) is the temperature-dependent conductivity 
of the material and E is the applied electric held. Explicitly neglecting the spatial variation 
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of T, Eq. can be simplified by a linear expansion in temperatnre: 


nC— = a{T)E^-a{T-To} (2) 

Here, Tq is the ambient temperature and a is an energy transfer rate from the sample to 
the environment depending on the specihc dissipation mechanism and the sample geometry. 
It has to be emphasized that inhomogeneous temperature and current distributions are 
neglected throughout this approach. Rather this model describes the electrical heating of 
a material in presence of just a heat flow out of the sample. By assigning the effective 
specihc heat to distinct microscopic processes of electronic or phononic (optical as well as 
acoustic) nature, the energy deposition into the individual microscopic subsystems and the 
dominant scattering mechanism of charge carriers can be analyzed. Therefore, an extended 
version of the model proposed by Mori is schematically illustrated in Fig. [T] explicitly 
taking into consideration the separation of the lattice system into the acoustic and optical 
phonon part and allowing for an additional charge carrier excitation channel via coupling 
to optical phonons. First of all, the energy provided by the electric held is fed into the 
electronic system. Depending on the electron scattering rate by acoustic and optical phonons 
(cte = tte-ac + ^e-opt), this energy can be transferred to the lattice system. However, if ae 
is small, hot carrier generation or direct population of excited electron states constitute the 
preferred relaxation mechanisms of the energy supplied by the electric held. This would 
be indicated by a low electronic contribution nCe to the specihc heat measured in this 
case. In this case, the temperature T = Te is assumed to be a parameter describing the 
population of excited electronic states rather than resembling the crystal lattice temperature 
Tl- On the other hand, identihcation of distinct contributions to the lattice specihc heat 
nCi = nCac + nCopt to nC yields evidence of an efficient carrier scattering by phonons. 
The relative magnitude of the electron scattering by acoustic and optical phonons also 
distinguishes the preferred excitation of the respective phononic subsystem. Yet, optical 
phonons do not contribute signihcantly to the heat transport because of their low group 
velocity, i.e. the energy stored in the optical phonon system is available quite long for an 
excitation of additional charge carriers while the heat stored in acoustic phonons quickly 
dissipates out of the system. Therefore, the measured a refers to limiting process of the 
energy transport out of the system. 

By numerical integration of Eq. the transient conductivity a{T\t]) of the material 
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FIG. 1. Illustration of the advanced electrothermal model. Depending on the charge carrier 
scattering rates, the electrical energy can be transferred to distinct microscopic subsystems of 
the material. An inefficient electron lattice scattering may lead to a direct excitation of charge 
carriers or hot electrons by the electric field. Due to their low contribution to thermal conductivity, 
energy stored in the optical phonon system is available for excitation of additional charge carriers 
while energy in the acoustic phonon subsystem quickly dissipates out of the sample. Therefore, 
identification of the microscopic contributions to the effective specific heat, the latter governing 
the electrothermal model, allows for an understanding of the energy flow in the system in relation 
to the preferred charge carrier scattering mechanism. 


for the applied electric field can be determined if the temperatnre-dependent condnctivity 
is known. The parameters a and nC can be adjnsted to match the simnlations with the 
measnred nonlinear current-volt age characteristics. To account for the observed nonlinear 
conduction phenomena in a variety of organic conductors, nC usually is assumed to be small 
compared to the specific heat of the material and is identified as the electronic specific heat 
nCe®®. Under these conditions, the temperature-dependence of the threshold electric field 
Eth as well as of the threshold current can be remarkably well reproduced. Nevertheless, 
a detailed analysis of the temperature-dependent effective specific heat as well as of the 
possible excitation of distinct lattice degrees of freedom accounting for the low nC observed, 
has not been presented so far. To account for voltage drops at the load resistor in our 
simulations we adjusted Eq. [^to the following form 
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dT _ { i+Rl%s(t) } rumoi -a 

dt ~ p-V-cUT)-Rs{T) p-V-cUT)^ 

■v* *v" 

P-term K-term 

with the sample volume V, the densitj^ p = 1.61 g/cm^, the temperature-dependent 
molar specihc heat CmiT) and the molar mass rrimoi- U is the voltage applied to the series of 
sample and load resistor {Rs and Rl) and 5 is the absolute energy transfer rate in \W/K]. 
The hrst term on the right hand side of Eq. [^will be denoted by P-term and second one by 
K-term from now on. Later, we also will allow for the specihc heat to become temperature- 
dependent. 

By neglecting the K-term and assuming a constant heating power {Rs =const), it is possible 
to extract the specihc heat from the linear relation AT = (P-term) ■ At on short time-scales. 
This method has already been proven to work reasonably well in measuring the specihc heat 
of rod-like samples, e.g. platinum wire^. Yet, the temporal resolution has to be sufficient to 
neglect the K-term which often is problematic in samples with higher resistance. Therefore, 
we approximate the K-term for short time-scales and small temperature steps by assuming 
a constant average temperature diherence 9 = T — Tq =const being half the temperature 
increase the transient linear ht extends to. Furthermore, we assume a constant ehective 
heating power by approximating Rs(T) = Rs,eff = .^5(^0 -|- 6 ) revealing a simple analytic 
solution of Eq. 


/\T = I ^^'^rnoli^ + RL/Rs,eff) ^ _ m^oi-a-d \ . 

\ pVcm{To+e)Rs^,ff pVc^{Ta+e) J '' 

{B-U^-C}-t = Afft ( 4 ) 

From the voltage dependence of the slope A^ff oc f/^ the specihc heat of the system under 
study can be determined quite reliably as the K-Term only represents a voltage independent 
ohset. We found that utilizing the specihc heat determined by this method allows us to 
simulate the transient resistance of DCNQI2CU very accurately. 

III. EXPERIMENTAL 

DCNQI-dg was prepared according to a literature procedure. High-quality single crys¬ 
tals of (DCNQI-d6)2Cu were grown electrochemically as described by Kato et al.^^Sl. The 
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needle-like crystals with typical cross-section areas of (0.03-0.1)^ mm^ and lengths of (0.2- 
20) mm were attached to free-standing gold wires of 25 pm diameter by silver-paint. Care 
has been taken to cover the ends of the crystal completely with silver paint to ensure a 
uniform current injection into the sample. The temperature-dependent resistance was mea¬ 
sured for different crystals from the same batch in two- and four-probe geometry in a He-flow 
cryostat between 4 K and 300 K with a Keithley 236 Source Measurement Unit (SMU). The 
data on nonlinear conduction presented here was measured in two-probe geometry on the 
same sample with a serial load resistor being roughly two orders of magnitude smaller than 
the sample resistance at low helds. Pulsed current-voltage characteristics were recorded by 
applying voltage pulses of 5-40 ms duration separated by 0.5 s intervals and measuring the 
current. In case of the transient resistance measurements, voltage and current were recorded 
with a Tektronix TDS744 oscilloscope. The applied voltage was fed into the oscilloscope by 
a 100:1 passive probe with =100 Mil input impedance while at the load resistor usually 
a 10:1 passive probe with =10 MU was used. A sketch of our transient setup is shown 
in Fig. 1^. 


IV. EXPERIMENTAL RESULTS 

A. Nonlinear Current-Voltage Characteristics 

The obtained crystals always revealed room temperature conductivities of 600-800 S/cm 
comparable to values in literatur d^ * ^^ ! The conductivity increased according to the typical 
a oc dependence down to about Tc = 85 K where its drop by hve orders of magnitude 
indicated the first-order Charge-Density-Wave (CDW) transition. Upon heating the transi¬ 
tion occurred with a small hysteresis at Tc = 89 K typical for hrst-order phase transitions. 
As depicted in Fig. i). the two-probe resistance measurement only reveals a jump by two 
orders of magnitude in resistance, i.e. the contact resistance can be neglected below the 
phase transition but dominates the conduction behavior above Tc. Above the phase transi¬ 
tion, the the contact resistance shows an activation energy of about 10 meV (fit not shown). 
Below Tc the resistivity is thermally activated with a corresponding energy of about AEa = 
(39±5) meV in agreement with values obtained from four-probe measurements on various 
crystals of the same batch. The interpolated two-probe resistance from the heating cycle of 
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Electric Field [V/cm] Time [10'^s] 


FIG. 2. (a) Setup used for the transient resistance experiments, (b) Two-probe sample resistance 
during the cooling and the heating cycle plotted logarithmically against the inverse temperature. 
The linear fit yields an activation energy of /\E k, 39 meV below the phase transition. The inset 
shows a microscope image of the sample, (c) Measured (dots) and simulated (solid lines) nonlinear 
current-voltage characteristics of (DCNQI-d 6 ) 2 Cu between 35 K and 75 K in steps of 5 K. (d) 
Simulated temperature transients after applying a voltage pulse of 87 V to the crystal at Tq = 40 
K. The transient P- and K-terms (Eq. contributing to the temperature rise dT/dt are illustrated 
separately. 

the sample with cross-section area of (0.03)^ mm^ and length of (0.275) mm was used for 
all subsequent simulations. 

Fig. illustrates the nonlinear current-volt age characteristics obtained from measurements 
with single voltage pulses of 40 ms length and 0.2-100 V magnitude at temperatures be¬ 
tween 35 K and 75 K together with the numerical simulation according to Eq. with Cm = 
55 J/(mole K) and a = 9T0“® W/K kept constant. The simulation is able to accurately 
reproduce the current-voltage characteristics and the temperature-dependence of threshold 





































field Eth and current jth- The reason for the reliability of the simulation can be understood 
by Fig. 1^ where the transient evolution of the temperature obtained from the simulation 
for a voltage pulse of 87 V at Tq = 40 K is shown. Here, already after 1.1 ms the steady 
state {dT/dt = 0) is reached which will be the case for most of the data points in the 
current-voltage characteristics. As the steady state solution of Eq. |^is independent of Cm, 
the final temperature and thus the resistance is determined by the choice of a as long as 
the value for is chosen not too high. Fig. also displays the individual contributions 
of the P- and the K-terms to the transient temperature rise dT / dt. The P-term governs the 
short-time behavior until the K-term reacts on the temperature rise leading to its steady 
state solution. It has to be noted that both terms still depend on Cm in our description, 
i.e. until the steady state is reached, the dynamic solution depends on Cm and 5. The 
applied value of Cm = 55 J/(mole K) corresponds to the specihc heat of (DCNQI-d 6 ) 2 Cu at 
30 hP. In order to gain detailed insights into the dependence of the resistive switching on 
microscopic contributions to the specihc heat, shorter pulses are required or, even better, 
the complete transient evolution of the sample resistance has to be analyzed which we will 
address in the following section. 


B. Dynamic Response 


The measured dynamic response of the sample resistance at Tq = 40 K to voltage pulses 
of t/ = 56 V and 17 = 97 V are exemplary illustrated in Fig. together with simulations 
assuming hxed values of and a. After a time-delay on the order of milliseconds, the 
resistance of the sample rapidly changes from a high-resistive to a low-resistive state. At 
high voltages a kink can be seen in the measurement data as well as in the simulations. It 
is caused by the huge drop of the resistivity at the phase transition and indicates when the 
transition temperature is reached, i.e. the charge order is melted. Adjusting the steady- 
state values of measurement and simulation after switching we determined a = (9 ± 1) ■ 10“® 


W/K in agreement with the value obtained in Sec. IV A The transient behavior is only fairly 
reproduced by the simulation, especially when the increase in sample temperature becomes 
large, i.e. on intermediate time-scales. This stems from an inaccurate determination of 
Cm = (55 ±25) J/(mole K) and its temperature dependence, the latter not being considered 
in the simulations. However, from the magnitude we can infer that the nonlinear conduction 
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in (DCNQI-d 6 ) 2 Cu is thermally driven, thus making it reasonable to treat the temperature 
parameter in the electrothermal model as the sample temperature. 



FIG. 3. (a) Measured (solid line) and simulated (dotted line) dynamic resistance response to 

voltage pulses of 56 V and 97 V and sample temperature (dashed line) calculated from the measured 
resistance at Tq = 40 K. At high voltages a kink is present in the measurement as well as in the 
simulations indicating when the phase transition temperature is reached, (b) Zoom-in on the 
initial temperature rise together with linear fits for different voltages. The slope Ag// of the 
linear fit reveals a (c) [/^-dependence from which the (d) specific heat can be determined for each 
temperature measured. In addition to the literature data by MatsuP^ and the expected acoustic 
phonon contribution, the fit of the effective Einstein model to our data is displayed. In analogy to 
the approach reported Matsui et al. we also added the latent heat of the phase transition (dashed 
line) to the Einstein model. 

In order to estimate the specihc heat and its temperature dependence more accurately 
and directly from our measurement data, we adopted the following evaluation procedure 
known from the transient electrothermal method^. First of all, we converted the transient 
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resistance into a transient sample temperatnre with the help of the temperatnre-dependent 
resistance cnrve measnred and displayed in Fig. |^. The resnlting cnrves are also illnstrated 
in Fig. 1^. Zooming-in the initial temperature rise we found a linear relationship of the 
sample temperature on short time-scales (see Fig. [^) as anticipated by Eq. The shaded 
area in Fig. corresponds to the overall temperature increase 26 in the htting regime. The 
expected tf^-dependence of the slope Aeff in the linear regime is depicted in Fig. |^. With 
Eq. I^we can determine the molar specihc heat from the slope B = by: 

“ pVBRs,,„ {1 + Rl/Rs,.,!? 

Here, Rs,eff = Rs(To + 6) represents the average resistance of the sample during the 
initial temperature rise 20 linearly htted with an error oi an = 1/2 ■ (i? 5 [To -|- 26] — i? 5 [To]). 
According to Eq. we additionally estimated a = (7 ± 4) • 10“^ W/K from the intercept 
C of the linear hts of Ae// for various temperatures. Given the low sensitivity of the initial 
temperature rise on 5, this value is in good agreement with the previously determined values 
of 5 = 9 ■ lO-*^ W/K. 

Including the error of Rs,eff and the slope B, the molar specihc heat Cm determined at 
different temperatures is presented in Fig. |^. The qualitative temperature dependence 
and the magnitude of agree quite well with specihc heat data of high accuracy published 
by MatsuP^. Nonetheless, within the error our measurement systematically undermines the 
published data by a rather constant ohset of about 40 J/(mole K). Even an inaccuracy in the 
optical determination of the sample volume cannot account for this deviation. In general, the 
specihc heat can be assumed to consist of contributions from acoustic and optical phonons 
as well as electrons: 


Cm Cac T Copt T Cel (6) 

The linear temperature dependence of our specihc heat data (not shown) with a slope of 
(3.0 ± 0.2) J/(mole K^) may be considered to originate from the electronic part Cei = 7 ■ T. 
Yet, the electronic specihc heat coefficient 7 = 0.025 J/(mole K^)ff^in the metallic state of 
DCNQGCu is two orders of magnitude smaller and therefore, should be negligible in the 
investigated insulating phase. Consequently, we disregard this contribution (and any other 
from low-energy charge-carrying excitation^2Ql)_ solids states with three translational and 
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three rotational degrees of freedom, the specihc heat contribution from acoustic phonons 
can be calculated by the Debye modeP^ 


a. = 18 i? I — 


’'Sd/T 




( 7 ) 


Jo (e^ - 1)^ 

R is the ideal gas constant and Qd the Debye temperature. Assuming a Debye tem¬ 
perature of ©ij = 82 this contribution to the specihc heat is also depicted in Fig. |^. 
According to the six degrees of freedom taken into account, the Debye contribution by acous¬ 
tic phonons approaches a high temperature limit of 6R ~ 50 J/(mole K). In spite of being 
on the same order of magnitude as our data, acoustic phonons only account for up to about 
44 J/(mole K) in the investigated temperature regime. Furthermore, the variation between 
35 K and 75 K is quite small due to the low Debye temperature of the material. Together 
with the observation that our specihc heat data is systematically lower compared to the 
published data by Matsui et al. this suggests, that the saturated contribution of low-energy 
acoustic phonons to the specihc heat does not allow for transient heating of our sample by 
inelastic scattering of charge carriers. Additional degrees of freedom are available by optical 
lattice phonons (external modes) and intramolecular vibrations (internal modes). Taking 
into account N = 22 atoms per molecule, there are 377 — 6 = 60 internal modes and 6Z —6 = 
30 external modes, Z = 6 denoting the number of building blocks constituting the unit cell 
of (DCNQI-d 6 ) 2 Cu single crystals. We estimated the contribution from optical phonons to 
the specihc heat in (DCNQI-d 6 ) 2 Cu by subtracting the Debye heat from the Matsui data 
and found a good agreement with our measurement data. Due to the low dispersion of both 
internal and external modes, their contribution to the specihc heat in molecular solids is 
reasonably well described by a series of harmonic oscillator^. To estimate the mean energy 
of the contributing optical phonon modes, we approximate our data by an ehective Einstein 
modeP^ 


c„, = N^R (^j _ 1)2 (8) 

where Ne is the number of oscillators with an ehective frequency ue = kB^E/h, and 
©F being the Einstein temperature. We obtained the best results by neglecting the elec¬ 
tronic {Cel) and acoustic phonon {Cac) contribution to the specihc heat and only htting 
the ehective Einstein model to our data, i.e. only taking into account internal and external 
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FIG. 4. Comparison of measured (solid lines) and simulated dynamic resistance response to (a) 
different voltage pulses at Tq = 40iF and (b) at different temperatures utilizing a constant specific 
heat of Cm = 55 J/(mole K) (dotted lines) and the effective Einstein model of the specific heat 


(dashed lines) as presented in Sec IV B 


optical phonon modes. The ht yields = 30 ± 3 and an effective Einstein temperature of 
Qe = (211 ± 9) K corresponding to an average energy of the contributing phonon modes 
of Eph = 19 meV. This mean energy correlates very well with the weighted optical phonon 
energies between 7-24 meV observed by Raman spectroscopy in DCNQRCrP^. Although the 
effective number of oscillators equals the number of possible external modes, it is most likely 
that internal and mixed modes contribute as welP^. Hence, we conclude that intra- and 
intermolecular optical modes govern the temperature dependence of the effective specihc 
heat in the investigated regime. To accurately simulate the complete dynamic resistance 
response we used the effective Einstein model of the specihc heat. The resulting transient 
simulations will be presented in the following. 


C. Simulation Dynamic Resistive Switching 


Because of the temperature raise by several ten Kelvin within the samples in response 
to the large voltage pulses applied we now explicitly take into account the temperature- 
dependence of Cm in the simulation (Eq. by means of our Einstein model described in 
Sec. IV B We also added a latent heat contribution (dashed line in Fig. |^) at the phase 
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transition to onr model, bnt fonnd its effect on the transients negligible in the temperatnre 
range nnder stndy as it only contribntes at large temperatnre differences between sample and 
environment, i.e. where the K-term already dominates Eq. i Fig. 1^ shows the transient 
resistance at Tq = 40 K in response to voltage pnlses of different magnitnde together with 
simnlations, both assnming a constant specihc heat of Cm = 55 J/ (mole K) (dotted line) and 
employing the temperatnre-dependent Einstein Model (dashed line). The transients are far 
better reprodnced inclnding explicitly the temperatnre-dependence of Cm in the simnlations. 
In particnlar, implementing the effective Einstein Model we can reprodnce the transients at 
varions temperatnres to a mnch better extent compared to the case of a constant specihc heat 
(Fig. 153 ). The small dip in the resistance curves after switching for temperatures above 50 
K is due to the current limiter of the pulsed voltage source which was used to prevent sample 
destruction. The overall good agreement between simulations and measurements underlines 
the accuracy of the specihc heat determined according to the procedure described in Sec. 

lYB 

V. DISCUSSION 

We ascribe the deviation of our data from literature to the respective measurement 
principle utilized. First of all, instead of steady-state conditioned heat how into and out of 
the sampl^i^, we measure the transient non-equilibrium internal heating of a sample caused 
by current how at intermediate to large electric helds. While we expect the temperature 
distribution inside the sample to be of minor relevance according to the short time scales 
being used to determine Cm, a general inhomogeneous current or held distribution inside 
the crystal might lead to an overestimation of the active charge transport volume in our 
analysis. The possible formation of an inhomogeneous high conduction state was already 
discussed by Mori and coworker^. This interpretation is supported by hicker noise studies 
revealing a pronounced noise level enhancement in organic charge-transfer salts which was 
also ascribed to the presence of an inhomogeneous current ho W^^ * ^^l Additionally, we were 
restricted to two-probe measurements due to the high helds necessary to induce the resistive 
switching. The non-planar silver paint contacts therefore might facilitate an inhomogeneous 
held distribution accompanied by local huctuations of the charge carrier density. 
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Nevertheless, being able to accurately model the specihc heat data by an effective Ein¬ 
stein model solely based on optical phonon modes, we rather explain the experimental 
hndings by extending the electrothermal model as shown in Fig. and allowing for an 
efficient stimulation of optical phonons by interacting with charge-carrying excitations. In 
contrast to acoustic phonons, the energy transferred to the optical phonon system remains 
within the sample as the contribution of these lattice excitations to the thermal conductiv¬ 
ity is small due to their low group velocity. This means that the population density of the 
optical phonon modes is shifted out of equilibrium and its temperature decouples from that 
of the acoustic phonon system. This excess energy facilitates the generation of additional 
charge carrier excitations leading to a multiplication effect. The mechanism of nonlinear 
conduction originating from current induced non-equilibrium optical phonon population has 
already been demonstrated in quasi-lD single-walled carbon nanotubed^^l as well as in 2D 
graphend^. 

The energy to excite additional charge carriers is estimated from the Boltzmann-like acti¬ 
vation of the DC conductivity with an energy of about Ea = (39 ± 5) meV slightly varying 
with temperature across the investigated regime (see Fig. [^). This energy is considerably 
smaller than the expected Peierls gap of Ecap ~ 100 meV and has been ascribed to an 
excitation of charge carriers from shallow states within the band gapP^. The activation en¬ 
ergy related to electronic transport is about twice as large as the observed effective phonon 
energy of Ep^ = 19 meV, for which reason two-phonon processes might be involved in 
the excitation of charge-carriers in the CDW state of DCNQI 2 CU. In general, an efficient 
interaction between charge-carrying excitations and optical phonons seems to be essential 
for the observed conductivity phenomena. 


Single excited charge carriers would need to gain the energy of Eph = 19 meV by 
the electric held to scatter with these optical phonons. Assuming an effective mass of 
rUe// = and an applied electric held of about E < 3600 V/cm, we estimate 

the necessary mobility of charge carriers to be /i > \J2Eph/meff/E = 1200 cm^/(Vs). 
This value is quite high but not unrealistic considering Hall mobilities on the order of 10® 
cm^/(Vs) which have been measured at 4 K for one-dimensional organic conductors such 
as (TMTSF) 2 PF^. Additionally, in the organic semiconductor naphthalene the saturation 
of charge carrier velocity at cryogenic temperatures and high electric helds has also been 
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ascribed to optical phonon generation via inelastic scatterin^Sl Alternatively, low-energy 
excitations of the electronic system might be constituted by soliton-like charged domain 
walls, which have been proposed to explain the low-frequency dielectric response in the 
CDW state of DCNQl 2 CiJ^. Soliton-like excitations are known to play a crucial role in 
the transport of commensurate CDW systems like polyacetylen^^ as well as of various 
organic charge-transfer salts, e.g. a-(BEDT-TTF) 2 l 2 P. It was noted that such solitons are 
expected to couple more efficiently to optical than to acoustic phonon modes by a factor 
of ^/a where a denotes the molecular lattice spacing along the chain direction and ^ is the 
spatial extension of the solitorpS. Thus, a large spatial extension of the soliton with respect 
to the lattice constant might lead to preferred interaction with optical phonons as indicated 
by our determined effective specihc heat. While we cannot completely resolve the funda¬ 
mental conduction mechanism from our data, an efficient interaction between the optical 
phonon modes and the charge-carrying excitations seems to be essential for the occurrence 
of nonlinear conduction effects in DCNQI 2 CU. 

For the sake of completeness, we also estimated the main loss mechanism of thermal energy 
to the environment. If the energy would be transferred from the crystal surface to the 
surroundings via convection or radiation, our determined a = 10“® W/K would yield a 
heat transfer coefficient of h = a/As = 3000 W/(m^K) with As = 3 ■ 10“'^ cm^ being 
the surface area of the sample. This heat transfer coefficient is much larger than that 
expected for convection [hconv = 6 — 30 W/(m^ K))l^. With the Stefan-Boltzmann con¬ 
stant as-, a maximum crystal surface temperature of Ti = 90 K and a minimum ambient 
temperature of Tq = 35 K, we can estimate the coefficient of radiative heat transfer by 
hrad = 0's ' {Ti + To){T/ -|- Tq) = 0.07 W/(m^ K) which, according to its magnitude, is 
also neglectecP^. The conductive heat transfer to the contacts is determined by the thermal 
conductivity k, of the material which we roughly estimate hj k < a ■ 1/{8 ■ A^s), I being 
the length and A^s the cross-section area of the sampld^. We obtain a value of k < 3.4 
W/(m K) which agrees reasonably well with the thermal conductivity of about 1.3 W/(m K) 
determined for alloyed DCNQI 2 Lio. 86 Cuo .14 at 50 K, bearing in mind the increase of thermal 
conductivity with increasing copper contenlpSl. Together with the assignment to acoustic 
phonons, the absence of appreciable variations of this thermal conductivity between 30 K 
and 150 K conhrms the interpretation of our temperature-independent a value. Hence, 
we identify the main energy loss mechanism in our samples to be the energy transferred 
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to the acoustic phonon system which is efficiently transported to the contacts via thermal 
conduction as depicted in our model (see Fig. [^. This energy exchange to acoustic phonons 
can occur via charge carrier-phonon scattering (ae-ac) or via scattering between optical and 
acoustic phonons {aopt-ac)- Both mechanisms seem to be inefficient at high electric helds 
supporting the occurrence of the measured nonlinear conduction effects. 


VI. CONCLUSION 

We experimentally studied the nonlinear conductivity in (DCNQI-d 6 ) 2 Cu single crystals 
by analyzing the dynamical resistance in response to large voltage pulses. In extension 
to simulations of the current-voltage characteristics at hxed pulse lengths, we demonstrate 
that applying the electrothermal description of low-dimensional conductor^ to the entire 
transient resistivity provides a more detailed insight into the energy exchange between the 
phononic and electronic subsystems of the material. Our simulations are able to accu¬ 
rately reproduce the nonlinear conduction behavior of our samples and reveal a distinct 
non-equilibrium excitation of optical phonon modes at a mean energy of about 19 meV. We 
infer a pronounced coupling between the electronic and optical phonon degrees of freedom 
at high electric helds in this material class. Due to their small contribution to the thermal 
conductivity, the energy transferred to the optical phonon system remains within the sample 
and can recover by additional single or collective charge carrying excitations leading to a 
thermally induced multiplication effect. Energy channeled into the acoustic phonon system 
is efficiently transferred to the contacts via heat conduction. As the involved effective optical 
phonon energy of around 19 meV is about half the activation energy of 39 meV observed 
in the DC conductivity, an excitation process including two phonons can be anticipated. 
Therefore, the nonlinear conduction in (DCNQI-d 6 ) 2 Cu crucially depends on the charge 
carrier interaction with optical phonon modes relative to the excitation of acoustic phonons 
via interaction of electrons or optical phonons with acoustic phonons. Its origin seems to 
be an interplay of population density of low-energy optical phonon states, small transport 
activation energy, large mobility at low temperatures and strong coupling between optical 
phonon modes and charge carriers at high helds in this material class. Hence, our results 
bridge the gap between the electrothermal model phenomenologically introduced by Mori 
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et al. and the microscopic origin of nonlinear transport effects in DCNQl 2 Cn. As snch, we 
consider the proposed mechanism also valid for nonlinear condnction phenomena in other 
low-dimensional molecnlar condnctors with different gronnd states. 
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